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We report on the numerical study of the relaxation rates of nuclear spins coupled through the 
hyperfine interaction to a two dimensional electron gas (2DEG) at magnetic fields corresponding to 
both fractional and integral Landau level (LL) fillings v. The Hamiltonians of up to 20 interacting 
electrons are diagonalized exactly in the spherical geometry, neglecting finite layer width, disorder, 
and LL mixing. The spectral functions t~ 1 (E) describing response of the 2DEG to the reversal of 
an embedded localized spin are calculated. In a (locally) incompressible v = 1 or | state, the finite 
Coulomb energy of short spin waves, together with the small nuclear Zeeman energy, prevent nuclear 
spin relaxation even in the limit of vanishing electron Zeeman energy (-Ez). However, we find that 
the nuclear spins can couple to the internal excitations of mobile finite-size skyrmions that appear 
in the 2DEG at sufficiently low Ez and at v slightly different from 1 or i. The experimentally 
observed dependence of nuclear spin relaxation rate on Ez and v is qualitatively explained in terms 
of the occurrence of skyrmions and antiskyrmions of various topological charge. 
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I. INTRODUCTION 

Two decades ago, transport experiments on the two- 
dimensional electron gas (2DEG) systems in a high mag- 
netic field B revealed rich new physics associated with the 
unique properties of their charge excitations, including 
the pair of most striking phenomena, the integralBo and 
fractionaHLl quantum Hall effects (IQHE and FQHE). 
Both effects are the manifestation of a finite gap opening 
in the charge excitation spectrum at a series of (integral 
or fractional) values of the Landau level (LL) filling fac- 
tor v = 1, 2, i, |, etc., and of the quasiparticle nature of 
the elementary charge excitations of this ft^fies of gapped 
(and thus incompressible) ground states a u 

Recent development of nuclear magnetic resonancec3 
(NMR) techniques allowed tbei^successful application to 
the quantum Hall systems,tirt2l and ultimately opened 
a new area of research associated with their spin de- 
gree of freedom. The quantum Hall systems with spin 
excitations are very attractive to both theory and ex- 
periment because of their fundamental aspects as well 
as for potential for applicatio ns ., n Th ey are liquids with 
unique (Laughlin) correlationsoEJIIj and unique excita- 
tions (iategrallyj-or fractionally charged quasiparticles 
withoutH or witlJl3 spin, skyrmionst3~E!3 and their_jCollec- 
tive excitations,EJ and chargeEj and spin wavesEj). On 
the other hand, the hyperfine coupling of the mobile elec- 
tron spin excitatipnS|-to the localized nuclear spins of the 
underlying atomsE-3c3 creates the possibility of control- 
ling the latter by inducing appropriate phase transitions 
in the 2DEG under variation of such experimentally ad- 
justed macroscopic parameters .as magnetic or electric 
fields, pressure, or temperature E3 If using nuclear spin 
states as a physical realization of quantum bits for storage 



of (and performing logical operations on) quantum infor- 
mation turned out to be possible, the coupled 2DEG- 
nuclei system should become a promising candidate, for 
the spin memory elements of a quantum computer .c3 

The most enlightening experiments on spin quan- 
tum Hall systems were those that offered NMR, evi- 
dence_for the occurrence of skyrmions in the IQHII3 and 
FQHlJ rtfSjpjC (confirmed by subsequent opticalo and 
transportE3E3 studies). Skyrmions (S^) and their con- 
jugates, antiskyrmions (S^), consist of K neutral spin 
waves bound to a particle in the empty reversecLspin (f) 
LL or to a hole in the filled (J.) LL, respectively.EHl In the 
IQH regime, the relevant particles are the reversed-spin 
electrons (cr) and the LL holes (h), and the skyrmions 
can be viewed as = e^.(eB,h)jc or = h{e#.ti)ic 
bound states, analogous to interband charged excitons 
Xf.B In theFQHpg ime, skyrmions consist of reversecU 
spin quasielectronsEJj (QEr) and Laughlin quasiholcscl 
(QH) bound to-brm S K = QE R (QE R QH) K and S+ = 
QH(QE R QH) X EJ The analogy between FQH and IQH 
skyrmions is most evident in the composite fermion (CF) 
pictureQ in which QEr and QH are represented by par- 
ticles and holes at the integral filling of their effective 
CF LL's. In both regimes, skyrmions are charged quasi- 
particles carrying large spin K. If the Zeeman energy 
Ei in a sample is sufficiently small for the isolated eR, 
h, QEr, and QH quasiparticles to become become un- 
stable towards creation and binding of a number (K) of 
spin waves to form skyrmions, the number of spin flips 
per particle added to or removed from a filled (electron 
or CF) LL is K. This quantity (K) sets the slope of 
the electron spin polarization (S z ), proportional to the 
Knight shift measured in the NMR experiments, as 
a function of the filling factor near v = 1, i etc. 
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Of the most recent ones, particularly-intriguing seems 
the NMR experiment of Kuzma et alE3 which revealed 
extremely long nuclear spin relaxation time in the rEQH 
regime, r < 0.5 s, exceeding times recorded earlierc£l by 
~ 10 3 . So different a relaxation time found in seem- 
ingly similar systems suggest that different microscopic 
mechanisms can be responsible for nuclear spin relax- 
ation, depending on experimentally variable conditions. 
Discussion of such mechanisms is the subject of this work. 

We report on detailed numerical studies of the hyper- 
fine interaction of the incompressible quantum Hall states 
at v = 1 and ^, as well as their spin excitations (spin 
waves), reversed-spin quasiparticles, and skyrmions), 
with the localized nuclear spins. The many-electron wave 
functions are obtained from exact-diagonalization calcu- 
lations carried out in the Haldane spherical geometry, 
neglecting disorder and excitations to higher orbital LL's 
or to higher quantum well subbands. The spectral func- 
tion t^ 1 (E) that describes response of the 2DEG to the 
reversal of an embedded localized spin and governs nu- 
clear spin relaxation time r for the particular microscopic 
2DEG-nucleus spin-flip process is calculated. 

We find that in a incompressible v = 1 or | state, 
the reversal of a nuclear spin creates a spin wave of a fi- 
nite wave vector k simply related to the area occupied by 
one electron, (k\) 2 w v (where A is the magnetic length). 
Since the spin wave dispersion Esw{k) begins at the elec- 
tronic Zeeman gap, -Esw(O) = Ez, the energy of a spin 
wave coupled to a nuclear spin exceeds Ez by a term 
Esw(k) — i?sw(0) which is of the order of the character- 
istic Coulomb energy E c — e 2 /A cx \[B, Since Ez and 
Ec are both much larger than the nuclear Zeeman gap, 
the energy conservation is expected to exclude creation 
(or annihilation) of spin waves as an efficient mechanism 
for nuclear spin relaxation. Consequently, very long re- 
laxation times r are expected for nuclear spins embedded 
in a (locally) incompressible IQH or FQH state, even if 
Ez could be made arbitrarily small (by means of appro- 
priate doping or application of pressure). 

At v slightly different from 1 or ^ skyrmions (or 
reversed-spin quasiparticles) appear in the incompress- 
ible quantum Hall liquid. The response function t^ 1 {E) 
is calculated for these objects and shown to have peaks 
corresponding to their "internal spin excitations" in 
which the skyrmion spin K and a ngul ar momentum L 
both change by one (Sk — > S/f±i).t3E3 It is also found 
that the oscillator strength r^ 1 of these transitions in- 
creases with increasing K . Since the energy gap for the 
internal skyrmion excitations is much smaller than Ez 
(and, in particular, it is equal to the nuclear Zeeman gap 
at the series of values of Ez), the skyrmion-nucleus spin- 
flip processes will be allowed by the energy conservation 
law, and provide efficient nuclear spin relaxation mecha- 
nism under experimental conditions. 

In both IQH and FQH regimes, our results imply criti- 
cal dependence of the nuclear spin relaxation rate on the 
presence of skyrmions in the 2DEG (dependent on v, Ez, 



well width, etc.), in agreement with experiments. How- 
ever, it is quite remarkable that (because skyrmions of 
the v = \ state occur only at much smaller values of Ez 
than skyrmions at v = 1) the energy required to create 
an electron spin excitation is much smaller in the IQH 
regime than in the FQH regime over a long range of Ez ■ 
This is opposite to the relation between charge excitation 
gaps, which scale as fuv c cx B for IQH states and a (much 
smaller) Ec oc \B for FQH states. 

At this stage of our study, we have limited ourselves to 
the calculation and analysis of a simple spectral function 
for the idealized many-electron states. This only allows 
for estimates of relative relaxation times r _1 for different 
microscopic processes, but not for their actual magni- 
tudes. The model also neglects the effects of (i) the finite 
width of the quasi-2DEG, the tilt of the magnetic field, 
the particular density profile g(z) in the direction nor- 
mal to the 2DEG plane, or.-the nuclear polarization pro- 
file £>n (•?) hi that direction;E£l (ii) the interaction-induced 
electron or hole scattering to higher LL's (i.e., LL mix- 
ing); (iii) disorder; or (iv) non-equilibrium processes. All 
these effects may become important in realistic experi- 
mental systems, leading to the reduction of the electronic 
interaction energy scale compared to the nuclear Zeeman 
energy due to the finite width of electronic wave func- 
tions (i) or screening (ii) , dependence of the nuclear spin 
relaxation rate on the correlation between g(z) and Qn(z) 
profiles (i), or localization of electronic excitations that 
relaxes the angular momentum conservation law in the 
spin-wave-nucleus or skyrmion-nucleus scattering pro- 
cesses (iii). We plan to study these and other possible 
effects in the future, and the motivation for the analysis 
of the ideal model used here is based on the fact that 
such model allows identification and classification of pos- 
sible elementary microscopic spin-flip processes, as well 
as the formulation of the involved selection rules which 
will be only modified to a various degree depending on 
the specific experimental conditions. 

II. MODEL 

A. Electron liquid 

As an extension of the earlier work on the spin p*<fi-i 
tations of the 2DEG in the quantum Hall regimejla€2l 
we study coupling of these excitations to the localized 
(e.g. nuclear) spins. The model to describe the 2DEG 
is that of Ref. |20|, except that it is now extended to in- 
clude the presence of a nucleus. In order to preserve 
the 2D symmetry of an infinite quantum well in a finite 
size calculation, the electrons are confined to a Haldane 
sphereQ of radius R. The magnetic field B normal to the 
surface is due to a Dirac monopolc in the center of the 
sphere. The monopole strength 2Q is defined in the units 
of flux quantum <j>Q = hc/e, so that 4irR 2 B = 2Qcj>o and 
A = R/^/Q is the magnetic length. The single-electron 
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states (monopole harmonics) are the eigenstates WB of 
magnitude (I) and projection (m) of angular momentum 
and of spin projection a, and they form 5-fold (g = 21+1) 
degenerate LL's labeled by n = I — Q. 

The cyclotron energy tku c oc B is assumed much larger 
than the Coulomb energy Eq — e 2 /A cx \jB. However, 
no assumption is made about the electron Zeeman en- 
ergy, and 7] = Ei/ Eq is a (small) free parameter of the 
model. As a result, only the a = —\ (J.) and +5 (t) 
states of the lowest (n = 0) LL need be included in the 
calculation, denoted simply by \ma). 

The many-electron Hamiltonian in the lowest LL is 

H = H c m lCT <4 2CT ' £ W' c «i4<7 (mim 2 |V| 

+ E c ™t c ™t^z. (1) 

where c\ na and c m(T are the electron creation and annihi- 
lation operators, the summation goes over all orbital and 
spin indices, and V is the Coulomb interaction potential. 
The iV-electron eigenstates are expanded in the basis of 
Slater determinants 

|mi<7i . . . m N a N ) = c\ ni(Ji . . . c\ nN(JN |vac) , (2) 

where |vac) is the vacuum state. Basis @ allows auto- 
matic resolution of two good many-body quantum num- 
bers, projections of spin (S z — ^2 Ui) and of angular mo- 
mentum (L z = ^rrii). However, the lengths of spin (S) 
and of angular momentum (L) are resolved numerically 
in the numerical diagonalization of each (S z , L z ) Hilbert 
subspace. The additional quantum number K — — S 
measures the number of reversed spins relative to the 
maximally polarized state. The many-electron states on 
a (finite) sphere converge to the states on an (infinite) 
plane in the Q = (R/X) 2 — ► 00 limit, only the spherical 
orbital numbers L and L z must be appropriately^ re- 
placed by the planar ones, the projections of total and 
center-of-mass angular momentum, M and Mqm- 

In the lowest energy states of the system described by 
Hamiltonian ([!]) near the integral or odd-denominator 
fractional filling of the lowest LL [v = 1 or |), a 
small number of spin waves or skyrmions move (to a 
good approximation, independently) in the appropriate 
incompressible quantum Hall "background" state. Be- 
ing charge neutral excitations, spin waves move along 
straight lines and carry linear momentum Tik. On a 
sphere, their linear orbits are closed into great circles, 
and the linear wave vector k takes on discrete values 
following from the quantization of angular momentum, 
L = kR = 0, 1, 2, ... . Skyrmions (Sk) on the other hand 
are charged, particle-like excitations that move along cir- 
cular cyclotron orbits similar to those of electrons. Their 
motion is therefore similarly quantized in both geome- 
tries, with the lowest LL of states having M = Mqm + K 
and Mqm = 0, ±1, ±2, . . . (on a plane) or L = Q — K 
and L z = L, L— 1, L — 2, ... (on a sphere), respectively. 

Both spin waves and skyrmions may become localized 
in the presence of disorder which has been ignored in 



this work. However, the dominant effects of such local- 
ization, at least in the weak disorder regime, can easily 
be predicted by analogy with the interband emission of 
neutral and charged excitons. For spin waves, the local- 
ization alters the relative occupation of different fc-states 
(in particular, it increases occupation of the k — state). 
However, we show later that spin waves do not couple ef- 
ficiently to nuclear spins regardless of k. For skyrmions, 
localization of their cyclotron orbits in the lower energy 
areas (without significant distortion of their wave func- 
tions) has two consequences: (i) Freezing of the positions 
of (few) skyrmions relative to (many) nuclei and thus 
variation of the electron spin polarization from one nu- 
cleus to another; this causes broadening of the Knight 
shift over the measured sample in the NMR experiment; 
derealization of skyrmions and restoration of uniform 
electron spin polarization at higher temperatures is called 
the "motional narrowing" of NMR lines, (ii) Variation 
of skyrmion energies due to confinement relaxes the en- 
ergy conservation law for the skyrmion-nucleus spin-flip 
processes, and thus broadens the minima of the nuclear 
spin relaxation time r as a function of Ez- Nevertheless, 
our main conclusions regarding the form of the continu- 
ous (due to spin waves) and discrete (due to skyrmions) 
parts of the 2DEG response function, and the role of the 
two types of spin excitations for nuclear spin relaxation 
remain valid independently of localization. 

B. Coupling to nuclear spin 

The coupling of the electron system (the "background" 
quantum Hall state and its spin excitations) to a single 
isolated localized nuclear spin will be described by the 
contact hyperfine interaction HamiltonianE2l 

F = Aj2^S k 6(r j -R k ), (3) 

where A is the coupling constant, and Sj and (r^ and 
Rfc) denote the spin (position) of the jth electron and kth 
nucleus, respectively. Moreover, the distance between 
nearest nuclei will be assumed sufficiently large to justify 
neglecting their direct dipolar interaction, and exclude 
any multi-nucleus phenomena. 

Due to the translational/rotational invariance of the 
2DEG, the position of the nucleus can be conveniently 
chosen at the north pole of the sphere, where all electron 
wave functions of the lowest LL vanish, except for \l "f) 
and |Z I). Then, ignoring the overall coupling constant 
(independent of the system size, R or Q) 7 the transverse 
part of F describing the spin-flip processes with AS Z = 1 
and projected onto the lowest LL simplifies to 

F = c Jf n c ll . (4) 

Clearly, the reversal of the localized (nuclear) spin is ac- 
companied by a reversal of an electronic spin spread over 
a (cyclotron) orbit of finite radius ~ A. 
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According to the Fermi golden rule, the oscillator 
strength for the transition between a given pair of 
initial and final electronic eigenstates, \i) and | f), is pro- 
portional to the square of the matrix element of F, 

r-/ = \(f\F\z)\\ (5) 

and, accordingly, the spectral function for a given initial 
(ground) state is 

r-\E) = rr\E) = £ T r/5[E (Ek - E f )}, (6) 
/ 

where E\ and Et are the energies of the initial and final 
states, respectively. Also, from Eq. (||), F couples the 
electron states with equal L z and S z different by one. 

Note that t^ 1 {E) is only defined up to a coupling con- 
stant that we are unable to calculate. Therefore, when 
comparing our results with experiment, only the rela- 
tive intensities of different microscopic spin-flip processes 
(e.g., processes that involve spin waves with different fc, 
skyrmions with different K, or skyrmions as opposed to 
spin waves) are meaningful, but we are unable to estimate 
the actual magnitude of the nuclear spin relaxation rates. 
However, only these relative intensities are a universal 
property of the electron quantum Hall system, virtually 
independent of many experimentally variable parameters 
(electron densitynauclear spin polarization profile across 
the 2DEG planeO, etc.). While spectral function (§) 
captures the essential physics of the interaction between 
an (ideal) Laughlin liquid and localized spins, the details 
of both the liquid and the spins must be included in a 
realistic calculation of the relaxation rates. 

Due to the simple form of the operator F in our ba- 
sis (||), the oscillator strengths can be easily evaluated 
for any known pair of |i) and |/) eigenstates. Therefore, 
the most difficult computational task is the accurate cal- 
culation of the many-electron eigenfunctions. Moreover, 
because of the breaking of both spatial and spin symme- 
try by the operator F, entire multiplets with different L z 
and S z must be calculated for each L and S. Hence, the 
proper identification of the relevant many-electron eigen- 
states that: (i) describe a planar system in the R/X — > oo 
limit, (ii) have significant r _1 , becomes essential. 

III. INTEGRAL QUANTUM HALL REGIME 

A. Spin waves 

We begin with the integral quantum Hall regime and 
the filling factor of precisely v = 1. We evaluate numer- 
ically the oscillator strengths r" 1 for all possible transi- 
tions induced by the operator F from the initial nondc- 
generate incompressible IQH ground state with L = 0, 
S z = —\N, and K = 0. From the commutator [F, S 2 } it 
can easily be shown that the spin-flip transition defined 
by F couples the K — initial state \i) to two different 



1 3 4 




2 4 6 8 10 1 2 0.0 0.2 0.4 0.6 0.8 1.0 
L E (e 2 IX) 



FIG. 1. (a) The spin wave energy spectrum (energy E 
vs. angular momentum L and wave vector k) at u = 1 calcu- 
lated for N = 20 electrons on Haldane sphere. The bubble 
diameters give the oscillator strength r _1 for the spin wave 
emission coupled to the nuclear spin reversal, (b) The re- 
sponse function of a planar v = 1 state to a nuclear spin 
reversal (oscillator strength r _1 calculated on a sphere for 
different N < 25 and normalized by \fN vs. energy E). (c) 
The response function r~ shown as a function of wave vec- 
tor k. Zeeman energy Ez is excluded and A is the magnetic 
length. 

subspaces, with K = or 1. However, the length of the 
total projection of F \i), the vector obtained by acting by 
F on the initial IQH ground state, onto the K = sub- 
space equals TV -1 in a finite iV-electron initial state, and 
thus it disappears in the N — > oo limit. Therefore, not 
surprisingly, the only spin excitations coupled to an infi- 
nite (planar) v = 1 state by F are those with K = 1. Of 
these, the only ones with significant oscillator strength 
t -1 turn out to be the spin wave states, which at the 
same time are the lowest energy excitations at v = 1. 

The numerical results for the spin-flip transitions cor- 
responding to creation (annihilation) of spin waves in a 
finite v = l state of N = 20 electrons at 2Q = N- 1 = 19 
are shown in Fig. |l](a) . In this and all other spectra in the 
paper E stands for the energy difference between the final 
and initial states and is given in the units of Eq — e 2 /A. 
The Zeeman energy Ez is not included. The horizontal 
axis shows the total angular momentum L (for the spin 
waves, L = kR), and the oscillator strength r _1 is pro- 
portional to the diameter of a bubble around each energy 
level marked with a cross. 

As could be predicted from Eq. (^), the reversal of an 
electronic spin induced by F occurs over an area corre- 
sponding to a cyclotron orbit. This sets a characteristic 
length scale £ ~ A for the efficient spin-flip process, and 
indeed in Fig. |l|(a) t _1 has a maximum at a finite L, 
while it vanishes in the limits of both small and large L. 

To determine the spectral function of a v = 1 state 
of an infinite 2DEG we have compared data obtained 
for different electron numbers, N < 25, and plotted the 
results together in Fig. |l|(b). The oscillator strengths 
for discrete values of E are multiplied by the factor 
^/N oc y/Q oc R/X, which comes from normalization of 
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the wave function of the extended spin wave over the 
entire sphere. All data points lie nicely on one regular 
curve that describes the spin wave creation/annihilation 
in both finite (spherical) and infinite (planar) systems. 
As it is the characteristic wave vector k (through the 
characteristic length scale £) rather than the energy E 
that determines the position of the maximum of t -1 , 
in the inset (c) we replot it as a function of k (only set- 
ting r _1 (0) to zero as appropriate for an infinite system). 
t -1 (k) is a more universal characteristic of the v = 1 IQH 
state than t~ 1 (E) in a sense that the spin wave disper- 
sion Egyjlk) , in an ideal 2D system derived by Kallin and 
HalperinJH3 in experiment may also depend on additional 
characteristics of actual 2DEG (e.g., the well width). As 
expected, T" 1 (fc) has a maximum at k ~ A -1 , which de- 
fines the spin-flip length scale of £ ~ A cx vl. 

Based on Fig. |l|, we make the following observations: 
(i) The incompressible v = 1 liquid responds to the rever- 
sal of a localized spin by emission of a spin wave whose 
kinetic energy, i?g W (fc), increases as a function of k; (ii) 
The response function t^ 1 (E) vanishes in both k = 
and oo limits, and it reaches maximum at E correspond- 
ing to k ~ A -1 ; (iii) When electron Zeeman energy is 
added, the energy of a k ~ A -1 spin wave that can cou- 
ple to a localized spin reversal is a sum of two terms, Ez 
and .Esw(A~Li-); and (iv) Since i?sw(A _1 ) w jEc (in an 
ideal systenS) is typically much larger than the nuclear 
Zeeman energy, the energy conservation prevents efficient 
relaxation through emission of spin waves in a (locally) 
incompressible liquid regardless of the value of Ez ■ 

This weak coupling of the v = 1 state to the nuclear 
spins obtained above agrees well with long nuclear spin 
relaxation tiroes measured at this filling factor at low 
temperatures.E3 It is evident from the experiments show- 
ing rapid increase of the relaxation rate when either v 
is moved awas^-from 1 or temperature is elevated (from 
2.1 to 4.2 K)li3 that charged excitations provide more 
efficient mechanism for nuclear spin relaxation than the 
spin waves. Let us then move on to the analysis of the 
spin-flip processes in the presence of such excitations, 
reversed-spin quasiparticles and skyrmions. 



B. Skyrmions 

It is well-known that an extra particle (a reversed-spin 
electron or a spin hole) added to the v = 1 ferromagnetic 
ground state induces and binds spin waves (whDse num- 
ber K depends on Ez/Eq) to form a skyrmionE£l (Sk), a 
particle-like charged excitations carrying (possibly large) 
spin K. In Fig. y(a) we replot the (anti)skyrmion energy 
spectrum calculated earlier in Ref. |2fJ for N — 12 elec- 
trons at 2Q = 12. Using the exact particle-hole symme- 
try within an isolated LL, this state can be mapped onto 
one in which an extra reversed-spin electron is added to 
the v — 1 state at the same value of 2Q — 12, and hence 
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FIG. 2. (a) The skyrmion energy spectrum (energy E vs. 
reversed spin number K) at v = 1 calculated for N = 12 
electrons on Haldane sphere. Open and full symbols corre- 
spond to the Zeeman energy Ez = and 0.02 e 2 /A, respec- 
tively (A is the magnetic length), (b) The energy gap A for 
skyrmion spin excitations (Sk — * Sk±i) at v = 1* as a func- 
tion of Ez, compared to the spin wave gap A = Ez at v = 1. 



cal axis gives the skyrmion energy E measured from the 
maximally polarized (K — 0) state which in this case cor- 
responds to one spin hole in the v = 1 state. On the hori- 
zontal axis wep&how the skyrmion spin (or size), K, which 
is also relatedEJ to its angular momentum, L = Q — K . 

The open symbols in Fig. |](a) mark the skyrmion en- 
ergy spectrum Eg(K), excluding the Zeeman energy. In 
an infinite system (and not only at v = 1, but also 
in the FQH regime discussed in the following section), 
it can be quite accurately approximated by Es(K) w 
-£[K/(K + l)] a . In an ideal 2D system, the bind- 
ing enecgy of an infinite skyrmion at v = 1 is known 
exactly ,Ea £ = \ Eq , and any choice of a ~ 1 cap- 

tures the most essential feature of Eg(K), which is that 
E S (K - 1) - E S (K) > E S {K) - E S (K + 1) for each K. 

Although our numerics yields a « 1.7 at !/ = 1, it 
is quite illuminating to solve the simplest case of a — 1 
(the equally simple arithmetics for a = 2 gives essentially 
identical answer; moreover, for v = | it actually seems 
that a w 1). When the Zeeman term KEz is added to the 
skyrmion energy as shown with full dots in Fig. |^(a) the 
ground state becomes a finite-size skyrmion with a cer- 
tain K (as marked with an open circle, for Ez — 0.02 Eq 
it turns out to be K = 3). Using our simple model 
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it will be denoted here as 



= 1= 



Again, the verti- 



we obtain that the transition between the Sk-i an d Sk 
ground states occurs at Ez = £/[K(K+l)) (in particular, 
the smallest skyrmion, Si, is stable below Ez = \£): 
and the excitation gap Ak from Sk to the lower of the 
Sk-i or Sk+i states [it is plotted in Fig. ||(b) for the 
numerical data of Fig. §(a)] reaches its maximum value 
of A K = £/[K(K + 1){K + 2)] at E z = £/[K(K + 2)]. 

As seen in Fig. ||(b), the most striking feature of the 
skyrmion energy spectrum in the presence of the (suffi- 
ciently small) Zeeman energy is that, in contrast to the 



5 



spin wave spectrum of the (locally) incompressible IQH 
liquid, the gap for spin excitations is greatly reduced 
compared to (and largely independent of) Ez- Indeed, 
it follows from our expression for the maximum of Ak 
that A K /E Z < (K+l)- 1 and A K /E Z < ^E z /£. More- 
over, the gap skyrmion Ak can be brought to resonance 
with an arbitrarily small nuclear Zeeman energy at the 
at entire series of Ez near the Sk—i Sk transitions. 

Being charged objects, skyrmions move along electron- 
like cyclotron orbits and repel one another through an 
effective short-range pseudepotential similar to that of 
electrons in the lowest LL.cj Such short-range repulsion, 
causes Laughlin correlations between the skyrmions ,Ej 
which therefore avoid high energy collisions with one an- 
other and behave as well-defined quasiparticles moving 
independently in the underlying IQH liquid. We know 
that in a similar, electron-valence-hole system, Laughlin 
correlations of interband charged excitons {X~) with one 
another as well as with surrounding electrons simplify the 
photoluminescence of a dilute [y < |)J£DEG to the re- 
combination of spatially isolated X~'s.Ej By analogy, we 
expect that the many-skyrmion effects can be excluded 
from the skyrmion-nucleus spin-flip scattering. 

To verify the possibility of efficient nuclear spin relax- 
ation through the interaction with skyrmions, we have 
evaluated the spectral functions r% (E) for the initial state 
\i) corresponding to a single skyrmion (or antiskyrmion) 
in a finite, iV-electron v = 1 system. Due to the finite size 
of the skyrmion, its ground state is a degenerate electron- 
like LL. On Haldane sphere, this LL is represented by the 
angular momentum multiplet at L = Q — K . Different 
values of L z = —L, —L+l, . . . , L label different cyclotron 
orbits, and the closest orbit to the north pole (the posi- 
tion of the nucleus) is that with L z — L for a skyrmion 
(S^) and L z = —L for an antiskyrmion {S^). Clearly, 
the value of Ti/ for a given i — > / transition involving 
a skyrmion depends on L z , which plays the role of an 
impact parameter of the skyrmion-nucleus collision. 

We have studied numerically a few systems with dif- 
ferent values of N and 2Q = N. For each N we have 
calculated the complete r, / spectra corresponding to the 
(anti)skyrmionic initial states \i) = Sq , , S% , . . . with 
different values of L z , and to all possible finite states |/). 
As an example, in Fig. || we show the results for N = 12 
and the initial states with K = 0, 1, and 2 (for K = 1 
data for two values of L z is shown) . The transition energy 
given on the horizontal axis is E = Ef — E^ (excluding 
Ez), and the units of r _1 on the vertical axis follow from 
Eq. (||). In the insets we display the corresponding en- 
ergy spectra in which the initial skyrmion states as well 
as all the final states can be identified for each transition. 

All the spectra are quite similar in that they can be 
decomposed into the quasi-continuous part at E > 
(here discrete because of size quantization) due to the 
response of the underlying IQH state discussed in the 
preceding subsection, and a single peak at E < due 
to a Sk ~^ Sk+i transition. Since we are mostly in- 
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FIG. 3. The spectral function of the hyperfine transition 
operator F (oscillator strength r _1 vs. energy E) at v = 1* 
calculated for iV = 12 electrons on Haldane sphere. Strong 
peaks at E < correspond to the "internal" skyrmion transi- 
tions Sk — * Sk+i indicated with arrows in the energy spec- 
tra shown in the insets (e-g). L z is the skyrmion angular 
momentum projection related to the impact parameter of the 
skyrmion-nucleus collision. Different frames correspond to 
different initial electron states: h (a), (b), and S% (d) 
close to the nucleus, and S± farther from the nucleus (c). 
Zeeman energy E-z is excluded and A is the magnetic length. 



terested in the processes that may conserve energy to 
allow efficient nuclear spin relaxation, let us neglect all 
E > transitions and concentrate on the Sk —> Sk+i 
one, whose energy including Ez can be made equal to 
the nuclear Zeeman gap. Clearly, it is only allowed for 
\L Z \ < Q — (K + 1) and its intensity quickly decreases 
when L z is increased from the minimum allowed value 
(such increase of L z corresponding to an increase of the 
nucleus-skyrmion average separation before the collision, 
i.e. of the impact parameter). Because we assume no lo- 
calization and thus allow that a skyrmion moves freely 
over the position of the nucleus (placed at the north pole 
by an arbitrary choice), it is more physical to consider 
t -1 summed over all allowed values of L z as a character- 
istic of this "internal skyrmion transition." 

In order to complete the analysis of the role of 
skyrmions in nuclear spin relaxation, to the facts that 
skyrmion-nucleus spin scattering has finite oscillator 
strength and that it can conserve energy, we ought to 
add the dependence on the skyrmion size, K. In Fig. |(b) 
we display the total (summed over all L z ) values of r _1 
for the Sk — > Sk+i transitions calculated for N = 12 
and plotted as a function of K. Surprisingly, the func- 
tion r^^if) has a maximum at a finite K. In order 



G 




N E z (e 2 IX) 



FIG. 4. (a) The oscillator strengths r _1 of the "internal" 
skyrmion transitions Sk — > <?K+i induced by a nuclear spin 
reversal, as a function of the inverse electron number, iV" 1 . 

(b) Same t , but calculated in a finite, N = 12 electron 
system and plotted as a function of the skyrmion spin K. 

(c) Same r _1 , but extrapolated to N — > oo and plotted as a 
function of the Zeeman energy Ez ■ The position of each peak 
is the value of Ez at which the transition energy A is zero. 
The solid line includes Gaussian broadening, (d) Same as (b) 
but for data extrapolated to N — > oo. 

to estimate t _1 in an infinite system, we have recalcu- 
lated r~ 1 (i ; r) for different N. The results are shown in 
Fig. ||(a) where we plot r _1 for K = 0, 1, 2, and 3 as 
a function of the inverse system size, iV -1 . Very regular 
dependence of t _1 on iV -1 for each K allows their accu- 
rate (quadratic) extrapolation to the iV -1 — > limit, as 
indicated with the solid lines. The result of the extrapola- 
tion is that in an infinite (planar) system, r _1 increases 
monotonously with increasing K, which indicates that 
nonmonotonic behavior in Fig. ^(b) is an artifact. 

The nearly linear increase of the extrapolated values of 
t _1 with increasing K shown in Fig. |](d) suggests that 
the total intensity t _1 (summed over all L z ) depends pre- 
dominantly on the skyrmion area. Because in experiment 
the differential cross-section for the skyrmion-nucleus 
collisions also depend on the skyrmion area, one can ex- 
pect that also the nuclear relaxation rate will increase as 
a function of K (at a constant number of skyrmions). In 
Fig. H|(c) we plot the r _1 peaks corresponding to sub- 
sequent Sk - > Sk+i transitions as a function of Ez at 
which the energy of this transition (E) vanishes (com- 
pare Fig. |2|(b) in which the closing of the gap is shown). 
Assuming that the skyrmion-nucleus spin scattering is 
a dominant nuclear spin relaxation process and that it 
is most efficient when E w 0, the curves obtained by 
broadening of the discrete peaks with Gaussians imitate 
the nuclear spin relaxation rate as a function of Ez- It 
is remarkable that when the value of Ez is lowered, the 
peaks with higher K are selected from the spectral func- 
tion, the separation between the neighboring peaks de- 
creases and their intensity increases. Let us stress that 
this expected behavior for v = 1^ differs qualitatively 
from what we predict at precisely v — 1 , where the relax- 
ation rate should monotonously increase with decreasing 



Ez and remain small even at Ez = 0. 

IV. FRACTIONAL QUANTUM HALL REGIME 

Recent NMR0 and optical experiments near v = \ 
revealed similar dependence of the electron spin polariza- 
tion on the magnetic field to that found earlier at v = 1 . 
In contrast to an earlier predictionj^j it now seems plausi- 
ble that the fast and weakly temperature-dependent nu- 
clear spin relaxation near v = -| is somehow related to 
the presence of skyrmions in the Laughlin liquid. 

Near v = (2p + 1) _1 (p >q1 is an integer) Laugh- 
lin correlationsfa allow mappingQ of the low-energy inter- 
acting electron states onto the noninteracting composite 
fermion (CF) states with an effective filling factor v* ps 1. 
The Chern-Simons transformation, in which 2p magnetic 
flux quanta are attached to each electron, results in the 
effective CF LL degeneracy of g* = g — 2p(N — 1). On a 
sphere,Q this replaces the electronic single particle angu- 
lar momentum I = Q rj \{^P + i)(N — 1) by an effective 
CF angular momentum I* = Q* ps ±(N — 1), where 2Q* 
denotes the effective CF monopole strength. 

There are two types of low energy charge-neutral ex- 
citations of Laughlin v l= i ground state, similar to the 
charge and spin waves£3 of the v = 1 state. This sim- 
ilarity lies at the heart of the CF picture|3 where these 
excitations correspond to promoting one CF from a com- 
pletely filled lowest (n — 0) spin-J, CF LL either to the 
first excited (n = 1) CF LL of the same spin (!) or to 
the same CF LL (n = 0) but with the reversed spin (|). 
Similarly to v = 1, charge and spin waves at v — ^ are 
composed of three types of elementary quasiparticles: a 
hole in the n = spin-| CF LL and the particles in the 
n = 1 spin-| and n = spin-| CF LL's, representing the 
Laughlin quasihole (QH) and quasielectron (QE) and the 
reversed-spin quasielectron (QEr), respectively. Each 
of these quasiparticles is characterized by such single- 
particle quantities as (fractional) electric charge, energy, 
or Landau degeneracy of the single-particle Hilbert space. 

Expecting similar behavior, we have carried out similar 
calculations for the v ~ \ filling as described in the pre- 
vious section for v w 1. Let us begin with the Laughlin 
incompressible FQH state at precisely v = h. We found 
that the reversal of a nuclear spin in this state creates a 
spin wave, in perfect analogy to what happened in the 
IQH regime. The main difference is that the spin wave 
at v = i consists of a QEr-QH pair whose interaction 
energy scale is about an order of magnitude smaller than 
it was for a e^-h pair at v = 1, predominantly due to the 
fractional QEr and QH charge (but also due to a larger 
size of the QEr and QH wave functions). In Fig. || we 
show the graphs for the FQH regime similar to those of 
Fig. |l|. The spectrum shown in Fig. ||(a) is for N = 6 
electrons at 2Q = 15. In the CF picture, this corresponds 
to six CF's filling exactly their lowest spin-| LL of degen- 
eracy g* = 2Q* + 1 = 6. Similarly to v = 1, additional 
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FIG. 5. Same as Fig. |l] but for v = j. More excited states 
appear in frame (a) in addition to the spin wave. 



weak transitions to higher states appear in a small sys- 
tem, but the spin wave will remain the dominant feature 
of the spectrum in the N — > oo limit. 

In Fig. ||(b) we overlay the spin wave spectra t^ 1 (E) 
obtained for different values of N. Similarly to v = 1, all 
data points (with r _1 multiplied by y/~N) fall on the same 
regular curve which, as expected, vanishes in both E = 
and E = oo limits, and reaches maximum at the energy 
E « 0.025 Eq, about an order of magnitude smaller than 
at v = 1. In Fig. |5|(c) we replot r _1 as a function of wave 
vector k = L/R. By analogy to v = 1, we expect that 
the length of spin wave most strongly coupled by F to 
a nuclear spin reversal corresponds to the smallest area 
containing one electron (and thus containing one unit of 
electron spin that must flip to compensate the nuclear 
spin). For the uniform v = | state, the average area per 
electron is 3 times larger than at v = 1 , yielding y3 times 
larger length scale £, and thus the maximum of r _1 (fc) is 
expected at k = ~ (V3A) -1 w 0.58A -1 . Indeed, this 
seems to be true of our Fig. [5](c). 

To summarize our results at precisely v = ^, the mech- 
anism of the coupling to the nuclear spins is very anal- 
ogous. The two major differences can be predicted from 
the simple arguments of 3 times reduced electron density 
and 3 times reduced charge of the involved quasiparticles. 
These differences are (i) about y3 times larger charac- 
teristic length scale of the response £ (wave length of the 
"active" spin wave) and (ii) about 3 2 smaller interaction 
energy i?sw(£ -1 ) of such active spin wave. While (i) is 
the reason for the reduction of t~ 1 of the corresponding 
transitions (compare the maxima of y/Nr^ 1 in Figs. |l| 
and ^]), (ii) should actually enhance nuclear relaxation 
due to spin waves (if Ez can be made disappear) as a 
result of the weaker violation of the energy conservation. 

Let us now turn to the spin-flip processes involving 
FQH skyrmions. Their energy spectra and gaps for the 
"internal" excitations are very similar to those at v = 1 
except for an overall reduction of the interaction en- 
ergy scalep-and breaking of the skyrmion-antiskyrmion 
symmctryHl In particular, graphs analogous to those in 
Fig. describe also the FQH skyrmions, only with about 
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FIG. 6. Same as Fig. ^|but for v = |. The initial states in 
different frames are: QH (a), QE (b), QEr (c), and (d). 



an order of magnitude smaller critical values of Ez ■ 

In Fig. |6| we display some of the t^ 1 (E) spectra cal- 
culated for N = 6 electrons at 2Q* = N and N + 2, 
corresponding to one skyrmion (or QE, or QEr) and an- 
tiskyrmion (or QH) in the Laughlin v = -| state, respec- 
tively. As in Fig. |^, in the insets we show the energy spec- 
tra in which the initial and final states for each transi- 
tion can be found. In perfect analogy to the IQH system, 
we identify the sequences of i^-induced transitions with 
E < that occur between the skyrmion or antiskyrmion 
states of different K: QE — > QE R — > — > S% —* ■ ■ ■ 
and QH -> S+ -> S£ -> . . .. 



V. COMPARISON OF IQH AND FQH REGIMES 

Let us finally compare the skyrmion and spin wave ex- 
citations and their possible coupling to the nuclear spins 
in the IQH and FQH regimes. It turns out that the fol- 
lowing behavior depicted in Fig. |t] that can be predicted 
from the simple arguments and existing analytical results 
alone is not far from the results of our exact calculations. 
The skyrmion energy spectrum is adequately reproduced, 
by Eq.' @, with E = \^J2E C for the IQH regime^ 
and an about 10 times smaller value for the FQH regime. 
This determines the dependence of the skyrmion size K 
and its gap A for spin excitations (Sk — > Sk±i) on Ez in 
both regimes. These plots of A(Ez) are marked as "IQH- 
Sky" and "FQH-Sky" in Fig. g On the other hand, it 
follows from the fact that A is the size of a cyclotron orbit 
that the spin waves must have k ~ A -1 to strongly colu 
pie to a localized spin reversal at v = 1. And knowingo 
the spin wave dispersion Esw{k) allows an estimate of 
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FIG. 7. The comparison of transition energies A of IQH 
and FQH systems corresponding to the spin wave emission 
and the internal skyrmion excitations, obtained from Eq. (Q) . 
Ez is the Zeeman energy and A is the magnetic length. Frame 
(b) shows a blow up of the shaded part of frame (a). 

the total (Zeeman plus Coulomb) energy gap for such 
"active" IQH spin waves, A ~ Ez + jEq- By reduc- 
ing the interaction energy by an order of magnitude one 
can also predict the spin wave gap in the FQH regime, 
A ~ Ez + jfiEc- These two plots of A(Ez) are marked 
as "IQH-SW" and "FQH-SW" in Fig. @ 

Since the energy conservation requires that A be equal 
or at least close to the nuclear Zeeman energy which is 
essentially zero, it is clear from Fig. [7] how the relative 
efficiencies of spin waves and skyrmions in both regimes 
depend on Ez ■ It is noteworthy that the (eR-/i or QEr- 
QH) interactions can have different effect on the spin gap 
A, depending on the presence of skyrmions in the system. 
Compared to a noninteracting system for which A = Ez, 
the interactions increase the spin gap A associated with 
the creation of spin waves, but decrease such a gap as- 
sociated with the skyrmion excitations. For spin waves, 
the enhancement of A is due to a decrease of e^-h or 
QEr-QH attraction at a finite wave vector k ~ ^Jv f\ 
(that can be interpreted as a spin wave kinetic energy). 
For skyrmions, the reduction of A is due to interaction 
induced level crossings and ground state transitions. 

It is clear from Fig. fjjthat the strong interactions in the 
IQH regime prevent the spin waves at v = 1 from com- 
ing into resonance with nuclear spins regardless of the 
value of Ez and practically eliminate them as an efficient 
nuclear spin relaxation mechanism at this filling. But 
at the same time, these interactions allow efficient relax- 
ation near v — 1 through the spin-flip nucleus-skyrmion 
collisions over a long range of Ez- On the other hand, 
the much weaker interactions in the FQH regime do not 
completely exclude nuclear relaxation by means of spin 
wave emission at Ez ~ 0, but they considerably shorten 
the range of Ez in which the skyrmions occur and can 
spin-flip collide with the nuclei. Therefore, a drop of the 
nuclear spin relaxation time r caused by the introduc- 
tion of charge excitations to the incompressible liquid (by 
varying density or magnetic field to move v away from 1 
or 3, increasing temperature, or inducing current) should 



be more pronounced in the IQH regime. This agrees with 
the experiments that typically show much longer relax- 
ation times at v = 1 than at v = ^. 

VI. CONCLUSION 

Using exact numerical diagonalization techniques we 
have studied possible relaxation mechanisms of nuclear 
spins coupled through the hyperfine interaction to the 
quantum Hall states of a 2DEG at filling factors near 
v = 1 and 4. By extrapolation of our finite-size results, 
we were able to determine the spectral function t^ 1 (E) 
describing response of an infinite (planar) 2DEG to the 
reversal of an embedded localized spin. We found that 
the spectral function can be decomposed into a contin- 
uous part describing transitions from the incompressible 
"background" state, and a discrete part which is due to 
the presence of additional charge excitations (skyrmions). 

The continuous part of the response function t^ 1 (E) 
describes emission of a spin wave, whose energy -Esw is 
a sum of the electronic Zeeman gap Ez and the kinetic 
energy dependent on the wave vector k. We found that, 
when expressed as a function of wave vector, r^ 1 vanishes 
in both k = and 00 limits, and that it has a maximum 
at a finite k m \~ X ^JD. Using the known spin wave dis- 
persion for v = 1 , we showed that the emission energy of 
an "active" spin wave that can couple to a nuclear spin 
reversal exceeds Ez by a kinetic term ~ \Eq ■ v 2 . This 
implies that even in the limit of vanishing Ez , the energy 
conservation law will prevent coupling of the (locally) in- 
compressible quantum Hall states to nuclear spins. This 
result agrees with long nuclear spin relaxation times ob- 
served in experiments in the absence of charged excita- 
tions (at precisely v = 1 or ^ and at low temperature). 

The situation changes dramatically when skyrmions 
are introduced into the 2DEG by either moving v away 
from 1 or |, increasing temperature, or applying volt- 
age to induce electric current. The reason for a dif- 
ferent behavior are the so-called "internal" spin excita- 
tions of skyrmions (in which their spin K changes by 
one) whose energy is much smaller than Ez and can be 
brought into resonance with the nearly vanishing nuclear 
Zeeman energy. Moreover, we have checked that the os- 
cillator strength t _1 of the skyrmion-nucleus collision 
corresponding to the K «-> (K + 1) transition is large 
and increases with increasing K. 

In both IQH and FQH regimes, our results imply crit- 
ical dependence of the nuclear spin relaxation rate on 
the presence of skyrmions in the 2DEG, in good agree- 
ment with the experiments. The contrast between the 
forbidden relaxation due to spin waves and the efficient 
relaxation due to skyrmions should be more pronounced 
at v = 1 than at v — i because of a larger (by an order 
of magnitude) interaction energy scale. 
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